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' We present the results of a lattice computation of the form factors for iv^l^vi decays near 

zero-recoil. These results will allow a determination of the CKM matrix element |Kib| when mea- 
surements of the differential decay rate become available. We also provide models for extrapolation 
' of the form factors and rate to the full recoil range. Our computation is performed in the quenched 

approximation to QCD on a 24'^ x 48 lattice at /3 = 6.2, using a non-perturbatively C'(a)-improved 
' action. The masses of all light valence quarks involved are extrapolated to their physical values. 

> . 
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Semileptonic decays of heavy-light mesons, particularly those containing a b quark, have attracted considerable 
interest. They play a crucial role in determining Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. With CLEO 
and now BaBar and Belle taking data, the prospects for an accurate measurement of the differential rate for B ttIvi 
decays are excellent^. Combined with a first-principles calculation of the relevant form factors, this measurement will 
allow a model-independent extraction of |Vub|, currently one of the least well-determined CKM parameters. 

■ The transition amplitude for exclusive semileptonic b ^ u decays factorises into leptonic and hadronic parts. The 
'— H hadronic matrix elements contain the non-perturbative, strong-interaction effects and are currently the largest source 

of theoretical uncertainty in the determination of |Kib| from these decays We present here a lattice QCD 

1^ , calculation in the quenched approximation of the matrix element relevant for tt^I'^vi decays. While this matrix 

■ element has already been obtained from the lattice at restricted values of |4| or for a range of with a "pion" 
^ ' composed of quarks with masses around that of the strange quark |^ , we determine it for the physical pion and a wider 

• range of values]^ We perform a fully C'(a)-improved calculation to reduce discretisation errors and implement light 

quark mass extrapolations as suggested in pO[ |. The q^-dependence of this matrix element has also been calculated 
^ using light-cone sum rules and a variety of quark models . 
- - ' The matrix element can be parameterised in terms of two form factors: 

(7r-(fc)|F^|i?"(p)> = f+{q^)(p + k - <zA„..)^ + /o(g')<z''A„,. (1) 

where q = p — k and A„j2 = {Mg — m^)/g^. The form factors /+ and /o are both real, dimcnsionlcss functions of the 
four-momentum transfer squared. In the limit of zero lepton mass the differential decay rate is given by 



0^ 
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'^CLEO already have a measurement of the total rate Q and have been able to obtain the differential decay rate for B ~* plvi 
decays in three recoil bins Q. They are currently attempting a similar binning for B — > irlvi decays 
^Preliminary studies of these decays have already been presented by two groups . 
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where A is a kinematic factor given by 

A(g2) = (Mi +ml- q^f - AMlml (3) 

and where q^, the invariant mass-squared of the lepton pair, takes values in the range from to {Mb — rnT^)"^. 

We work in the quenched approximation to QCD on a 24"^ x 48 lattice aX (3 — 6.2, corresponding to an inverse lattice 
spacing of = 2.bAt,\ GeV as determined from the p mass iQ. We use 216 gauge configurations generated from 
a combination of the over-relaxed and Cabibbo-Marinari [|l6| algorithms with periodic boundary conditions. 

The heavy-light mesons are composed of a propagating heavy quark with mass around that of the charm quark to 
keep discretisation errors in check, and a light antiquark with a mass around that of the strange quark to keep finite 
volume errors and the CPU time required to obtain the corresponding quark propagator under control. We then 
obtain results for t he p hysical Tr~l~^ui decay by extrapolating in heavy and light quark masses. We employ an 

SW fermion action jl^ that removes all 0{a) discretisation errors from the hadron spectrum We can also remove 
all 0{a) discretisation errors from matrix elements of local currents by an appropriate definition and renormalisation 
of these currents iflql. In the case of the vector current, we have 



= y^: + cy a- (a, + di)Ti^ (4) 



where l/f = Qlp.q and T^^^ = Qcrf^^q are the local lattice vector and tensor currents respectively and where di, {d*) 
is the forward (backward) lattice derivative. The renormalised current is 

K ^Zv(^l + ^-{amQ + amS^ + O{o?) , (5) 

where amQ ^ — l/2KQ.q — l/2Kcrit and both by and the matching coefficient Zy have been determined non- 
perturbatively p9[ |. The mixing coefficient cv is only known to one loop in perturbation theory [2C| ], but is small. 

To quantify residual discretisation errors, we consider the effective matching coefficient, Z^ , obtained from the 
forward matrix element of the current, — Qj^Q, between degenerate heavy-light mesons at rest. Since, in the 
continuum, this current is conserved we have 



(P(0) \Vo\P{0) ) = 2Mp = Zf(Pi6) \V,'\P{6) ) (6) 

where P is a heavy-light pseudoscalar meson. We calculate Zy^ for two values of the heavy quark mass. We 
compare these results with the matching coefficient, Zy(l -j-byomg), evaluated for these quark masses using the non- 
perturbative values of the coefficients Zy and by- These two measures will differ by terms oiO{a^). The comparison is 
shown in Table |. The agreement between the two procedures is excellent, suggesting that, at least for short-distance 
quantities, higher-order discretisation errors are small. However, we cannot exclude the possibility that discretisation 
errors for non-degenerate, non-forward matrix elements are larger than the above comparison suggests. Although the 
present calculation has been performed at one value of the lattice spacing, recent work [ pT| using the same improved 
action, suggests that for long-distance quantities such as masses the error in extrapolating to the continuum limit 
from (3 — 6.2 is around 6%. 

Four values of the hopping parameter, k, were chosen for the heavy quark {kq — 0.1200, 0.1233, 0.1266, 0.1299), 
three values for the active light quark (A) (ka — 0.1346, 0.1351, 0.1353), and two for the spectator light quark (S) 
{ks — 0.1346, 0.1351). We obtain the relevant matrix elements from heavy-to-light three-point correlation functions, 
divided by the appropriate factors extracted from two-point functions as in Refs. |^,^. We place the operator for 
the heavy-light pseudoscalar a,t t = 20, and the operator for the light-light pseudoscalar at t = 0. We find that 
contamination of the signal from other time orderings as well as from excited states is negligible. 

We use eight different combinations of p and k to determine the momentum dependence of the corresponding form 
factors: ^ 0, ^ 1, ^ v^, 1 ^ 0, 1 ^ 1, 1 ^ 1±, 1 ^ 1^ and 1 V2i_ in lattice units. There is no 
channel for /_|_. For each momentum channel, the form factors are obtained from a simultaneous fit to the matrix 
elements of the temporal and spatial components of the vector current. 

The chiral limit is reached at the value of the hopping parameter for which the pion mass vanishes. For the present 
simulation this value is Kcrit = 0.13582(1). The physical value of m^/mp is reached for k„ = 0.13578(1). In order 
to evaluate the form factor, /i, i = +, 0, at physical light quark masses, we must consider both the explicit mass 
dependence of fi and the indirect dependence arising from the change in q^: 
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f^ = Ml{^A,^s)^'^A,Ks). (7) 

We first interpolate, to a common set of values, the form factors corresponding to a given heavy quark and different 
light quark mass combinations. For the different heavy quarks, these sets are chosen such that the corresponding sets 
of heavy quark velocities are the same for all heavy quarks, for reasons that will become clear below. The interpolation 
is performed at fixed light quark mass using a form motivated by pole dominance models: 

^ ■ («' 

with the kinematical constraint, /+(0) — /o(0), imposed. Any consequent model dependence is mild as we use this 
ansatz only to interpolate the form factors in the range of for which we have data. We have also tried other ansatze 
including dipole/pole for f+{q^) and fo{q^) respectively, with and/or without the kinematical constraint. This is 
included in our estimate of systematic errors and is shown in Figure ^ 

These interpolated data points are then extrapolated in ka and Kg to k„ at fixed q^ , according to 

f{ns,KA) = a + (3(—~—\+-i(— + - —\ (9) 

The extrapolation is two-dimensional in the active and spectator light quark masses. 

In previous UKQCD collaboration analyses |^,^ a term linear in the pion mass was added to account for the 
indirect q^ contribution. Here we hold q^ fixed as the light quark masses are changed ||lo|] . This yields a more reliable 
extrapolation by separating the q^ and the explicit light quark mass dependences. 

Heavy quark symmetry (HQS) is used to model the form factors' dependence on heavy meson mass at fixed four- 
velocity, V. We thus work with the recoil variable 

^.,^M^±<^^ (10) 
2Mp ^ ' 

The scaling relations for /-|_ and /q at fixed v ■ k, as given by heavy quark effective theory are p^ , 

C(Mp, Mb)Mv ■ fc)A4*/' ^ ^.^ (^1 + A + _^ + . . (11) 

where the ellipsis denotes higher order terms in the heavy quark expansion and Si = {—1, 1} when i — {-f-,0}. The 
coefficient C is the logarithmic matching factor p5[|. 
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where /3o = H in quenched QCD and is the one loop running coupling with ALl. = 295 MeV. 

We choose a set of five values of • A; in a region accessible to our calculation for all of our values of heavy and light 
quark mass. For each initial heavy quark this yields a distinct set of q^ values. It is these sets of q^ values which we 
use in the light quark mass extrapolation described above. From that extrapolation, we obtain the form factors for 
all Mp at the chosen values oi v ■ k. We then extrapolate /+ and /q at fixed v ■ k to obtain the form-factors at the 
B-meson mass. Figure |] shows both a quadratic fit to all four data points and a linear fit to the three heaviest. The 
difference between these two procedures is used to estimate the systematic uncertainty in this extrapolation. Table || 
shows the heavy quark and meson masses used in this work. 

The resulting 7r~Z+i/; form factors and differential decay rate are shown as functions of q^ in Figures ^ and ^ 



and their values are summarised in Table [II . It is important to note that these results can be used for a determination 
of |Viib|, without model-dependent assumptions about the q^ dependence of form factors, once experiments measure the 
differential or partially integrated rate in the range of q^ values reached by our calculation. Future lattice calculations 
in full, unquenched QCD will permit a completely model-independent determination of |Vub|- 

The central values and statistical, bootstrap errors in the figures and table are obtained with the kinematically 
constrained pole/pole g^-interpolation function of Eq. (^, quadratic heavy-quark-mass extrapolations and with the 
lattice spacing set by rrip. Systematic errors are determined by considering the change in our results due to: different 
choices of g^-interpolation function, as discussed after Eq. (||); using tq |^6| to set the lattice spacing; and using a linear, 
instead of quadratic, heavy extrapolation for the heaviest three quarks. Using the KLM quark field normalisation |2^] 
produced negligible change in these results, suggesting that higher order mass-dependent discretisation effects are 
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small. We also assign a 6% error representing the possible effect of a continuum extrapolation as mentioned above. 
The quoted error is the quadratic sum of the maximum variation from all sources of systematics. 

To extrapolate these results to a larger range of g^, we make model assumptions |^ Pole dominance models suggest 
the following momentum dependence for the form factors, 

where i = +, 0, is an integer exponent and the kinematical constraint /+(0) = /o(0) has already been imposed. 
Combining this with the HQS scaling relations of Eq. ( JTl] ) implies n+ = uq + 1. Light-cone sum-rules scaling further 
suggests no = 1 ^. Another pole/dipole parameterisation for /o and /+, which accounts for the B* pole in /+ 
correctly, has been suggested by Becirevic and Kaidalov (BK) pG|: 



f+iq 



(l-gVA/|0(l-"'?VMi 



We fit both models to our form factors extrapolated to the B-meson mass. The results are shown in Figure || for 
our central value procedure. Though both parameterisations fit the data equally well, we favour the more physical 
BK description. 

These models can then be used to compute the total decay rate by integrating Eq. (|^) with respect to q^. This rate 
of course depends on the choice of model and is sensitive to uncertainties in the lattice results, since it is dominated 
by the low region, only reached by extrapolation. We find 

r/lKbP - {9tltl) ps-\ (15) 

The first error is statistical and the second is the systematic error estimated, as for the form factors, from the effects 
decribed above as well as the difference between the two models. Combining this result with CLEO's exclusive 
iKbl = (3.25 ± O.Utoil ± 0.55) x IQ-^ § and the B° lifetime, Tgo = 1.54 ± 0.03 ps |l|, we find a branching ratio 
r/Ftotai = (l.5lg;3tQ;3 ± 0.6) X 10~*, where the first error is the lattice statistical error, the second systematic and the 
third, the experimental errors on the mean combined in quadrature. This result is consistent with the measurement 
by the CLEO collaboration 0, r/Etotai = (1-8 ± 0.4 ± 0.3 ± 0.2) x 10"*, where the errors are statistical, systematic 
and from model dependence. 

In this letter we have reported a lattice computation of the form factors and differential decay rate for the decay 
n^l^vi. We have used a fully C'(a)-improved action to minimise discretisation errors. We are repeating the 
calculation at a different lattice spacing to probe the effect of improvement. It should also be remembered that 
the computation was performed in the quenched approximation; we are currently generating configurations with two 
flavours of light dynamical quarks and aim to quantify the quenching effect. 
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FIG. 1. Interpolations of f+{q^) and fo{q^) for the heaviest of the heavy quarks and the lightest of the light quarks. The 
support of the curves shows the range of which results in interpolation for all light quark mass combinations. The fits shown 
enforce the kinematic constraint /+(0) — /o(0). 




FIG. 2. Heavy-quark-mass extrapolation of f+{q^) at fixed v ■ k corresponding to q^ — 22.3 GeV^ at the B meson scale. 
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FIG. 3. Momentum dependence of the form factors. The data shows statistical errors only. 
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TABLE I. The effective matdiing coefficient for different values of tlie heavy quark mass and spectator hopping parameter. 
The current does not depend on the spectator quark. 

amg Zf, Ks = O.f 346 Zf, ks = O.f 35i Zv{l + hvamq) 

0.485 f.3f5t^ f.Sfelg f.332t^ 

0.268 imitl f.085^2 f.ogotg 



TABLE 11. Heavy quark and meson masses used in this work. The table shows the hare quark mass in lattice units, the 
renormalisation group invariant quark mass defined by — Zm{1 + hMamQ)mQ |33] and the heavy-light (k„) pseudoscalar 
meson mass. The scale is set by rrip and hm is obtained using the one-loop result of |2(| 

K amg (GeV) Mp (GeV) 

0.1200 0.485 1.52 2.035 ± 5 

0.1233 0.374 1.23 1.771 ± 5 

0.1266 0.268 1.02 1.483 ± 5 

0.1299 0.168 0.69 1.157 ±5 



TABLE IIL Form factors and differential decay rate as functions of (f . The central value comes from the constrained 
pole/pole interpolation procedure, the first error is statistical, the second is systematic as described in the text. 

{GeVf 16.7 18.1 19.5 20.9 22.3 

M^) 0-5711 O.eit^ljg 0-66^ 111 0.72tlt,l 0.79^1 

l/|Kbpdr/dg2 (ps-^GeV-^) 0.29^1° t'l 0.27t? t^l 0.25tl t'l 0.23tl t'l 0.20l^ t I 
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